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ABSTRACT 

We report the discovery of a large-scale coherent filamentary structure of Lya 
emitters in a redshift space at z = 3.1. We carried out spectroscopic observa- 
tions to map the three dimensional structure of the belt-like feature of the Lya 
emitters discovered by our previous narrow-band imaging observations centered 
on the protocluster at z = 3.1. The feature was found to consist of at least three 
physical filaments connecting with each other. The result is in qualitative agree- 
ment with the prediction of the 'biased' galaxy-formation theories that galaxies 
preferentially formed in large-scale filamentary or sheet-like mass overdensities in 
the early Universe. We also found that the two known giant Lya emission-line 
nebulae showing high star-formation activities are located near the intersection 
of these filaments, which presumably evolves into a massive cluster of galaxies in 
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the local Universe. This may suggest that massive galaxy formation occurs at 
the characteristic place in the surrounding large-scale structure at high redshift. 

Subject headings: cosmology: observations — cosmology: large-scale structure of 
the universe — galaxies: high-redshift — galaxies: formation 

1. INTRODUCTION 

Theories of structure formation predict that galaxy formation preferentially occurs along 
large-scale filamentary or sheet-like mass overdense regions in the early Universe and the 
intersections of such filaments or sheets evolve into dense clusters of galaxies at the later 
epoch (Governato ct al. 1998; Kauffmann et al. 1999; Cen & Ostriker 2000; Benson et al. 
2001). Recent deep observations of star-forming galaxies at high redshift, such as Lyman 
break galaxies (LBGs) or Lya emitters (LAEs), have revealed their inhomogeneous spatial 
distribution over several tens to a hundred Mpc (Steidel et al. 1998, 2000, hereafter S98, 
SOO; Adelberger et al. 1998; MoUer & Fynbo 2001, Shimasaku et al. 2003, 2004; Palunas 
et al. 2004; Francis et al. 2004; Ouchi et al. 2005). However, their volume coverage or 
redshift information are still limited. Thus, there is little direct observational evidence of 
this theoretical prediction. 

A large overdensity of LBGs and LAEs was discovered at 2; = 3.1 in the SSA22 and 
it was regarded as a protocluster of galaxies (S98 and SOO). We carried out wide- field and 
deep narrow-band imaging observations of the SSA22 proto-cluster and discovered a belt- 
like region of high surface density of LAEs with the length of more than 60 Mpc and the 
width of about 10 Mpc in comoving scale (Hayashino et al. 2004, hereafter H04). We could 
not distinguish, however, whether the belt-like feature on the sky is a physically coherent 
structure or just a chance projection of isolated clumps in space. 

There exist two giant extended Lya nebulae (Lya blobs, LABs) in this proto-cluster 
whose sizes are larger than 100 kpc (SOO). We also detected in our narrow-band images 
33 LABs with Lya isophotal area larger than 16 arcsec^, which corresponds to 900 kpc^ or 
d 30 kpc ai z — 3.1 (Matsuda et al. 2004, hereafter M04). It is, however, noted that 
M04's LABs are smaller than SOO's two giant LABs. The two giant LABs seem to be rare 
outstanding objects in the region. 

In this letter, we present the redshift distribution of LAEs in the belt-like feature. We 
use AB magnitudes and adopt a set of cosmological parameters, = 0.3, ^Ia = 0.7 and 
Ho = 70 km s-i Mpc-^ 
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2. OBSERVATIONS 

The targets of our spectroscopy were chosen from the 283 candidate LAEs detected in 
our narrow-band imaging observations centered on the SSA22 proto-cluster (RA=22'^17"^34^.0, 
Dec=+00°17'01", [J2000.0]) at z = 3.1 using the Prime Focus Camera (Suprime-Cam, 
Miyazaki et al. 2002) of the 8.2 m Subaru telescope. We briefly describe the selection 
criteria below. Details are given in H04. The 283 candidate LAEs satisfies the following 
criteria; (1) bright narrow-band magnitude at 4970A, NBA97 < 25.8 mag, and (2) large ob- 
served equivalent width, EWohs > 154 A (or L(LyQ;) > 10^^ erg s~^ at z = 3.1). In addition, 
we imposed one of the following two criteria; (3) red continuum color, B — Vc > 0.2 mag, 
for the objects with Vc brighter than 26.4 mag, or (4) larger equivalent width, EWqTos > 267 
A, for the objects with Vc fainter than 26.4 mag. Here, Vc represents the emission-line free 
V-band magnitude obtained after subtracting the narrow-band flux from the V-band flux. 
Note that LAEs at 2; ~ 3, which often have continuum spectra of fi, ~ should have 
B — Vc greater than 0.2 mag, since the Lyo; forest dims the continuum blueward of Lyo; line 
of galaxies (e.g. Madau 1995). 

The distribution of the 283 LAEs showed the belt-like feature of high surface density 
running from NW to SE on the sky (the bottom right panel of Figure 1). In order to examine 
the three-dimensional structure of this belt-like feature, we carried out the spectroscopic 
observations of the candidate LAEs using Faint Object Camera And Spectrograph (FOCAS, 
Kashikawa et al. 2002) on the Subaru telescope in the two first-half nights of 29 and 30 
October 2003 (UT). The diameter of the field of view of a FOCAS slit mask was 6'. We 
used 6 slit masks along the belt-like feature (see Figure 1). Out of the 283 LAEs, 122 are 
located in the area covered by the 6 masks. We prepared the slits for bright 84 of the 122 
selected in order of NB magnitudes. We used the 300B grism and the SY47 filter. The 
typical wavelength coverage was 4700-9400 A. The sht widths were 0".8 and the seeing was 
0".4-l".0. The spectral resolution was FWHM ~ 10 A at 5000 A, corresponding to a velocity 
resolution of ~600 km s"^. The pixel sampling was 1.4 A in wavelength (no binning) and 
0".3 in spatial direction (3 pixel binning). Exposure times were 0.7-1.5 hours per mask. We 
reduced the data using the IDL program developed by FOCAS team and IRAF. We used 
both Th-Ar lamp spectra and night sky lines in wavelength calibration. The rms of fitting 
error in wavelength calibration was smaller than 2 A. We used the standard star Feige 110 
in flux calibration. 
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3. RESULTS AND DISCUSSION 

Among the 84 targeted LAEs, 56 spectra show a single emission-Une with a signal-to- 
noise ratio (S/N) per resolution element (10 A) larger than 5. The typical S/N of the 56 
emission-lines is about 11. The most plausible interlopers of the emission-lines are [OII]A3727 
ai z — 0.325 — 0.346. In this case, the [OIII]AA4959,5007 should be detected, since the wave- 
length coverage of most spectra extends longer than 6700A. Most of star-forming galaxies at 
low redshifts have [OIII]A5007/[OII]A3727 ratio larger than 0.15 (e.g. Jansen et al. 2000). 
Since the upper limits for [OIII]A5007/[OII]A3727 ratios are smaller than 0.15 for our spectra, 
there is no evidence for contamination of [Oil] emission-line galaxies. Therefore, we identify 
the single emission-lines as Lya with high confidence. The mean redshift of the 56 identified 
LAEs is < 2; >= 3.091 and the redshift dispersion is Az = 0.015 {Av = 1100 km/s). We 
stacked the observed frame spectra of remaining 28 unidentified objects with emission-lines 
with S/N< 5. The stacked spectrum shows a significant double-peak emission-hne, whose 
profile is similar to the shape of the redshift histogram of the 56 identified LAEs. Accord- 
ingly, it is highly probable that a large fraction of the unidentified LAEs is also located in 
the same structure at 2; ~ 3.1. 

Redshifts of the LAEs are not the exact measure of the Hubble fiow or the spatial 
distribution due to their peculiar velocities. In fact, peculiar velocities are considered to be 
the dominant source of errors in the spatial distribution of LAEs^ . While the peculiar velocity 
dispersion of galaxies is 500 — 600 km s~^ in the local universe (Zehavi et al. 2002; Hawkins 
et al. 2003), it is expected to be smaller at high redshifts even in the over-dense regions 
(Hamana et al. 2001, 2003). Indeed, the predicted pair wise peculiar velocity dispersion of 
galaxies at 2; ~ 3 in cosmological simulations is 300 — 400 km s~^ (Zhao et al. 2002), which 
corresponds to a very small redshift dispersion of ~ 0.005. 

In Figure 1, we plot the resultant three dimensional distribution of the 56 LAEs, together 
with the projected distribution. We can see that the belt-like feature seen in the projected 
distribution consists of a filamentary structure running from (ARA[arcmin], ADec[arcmin], 
Redshift [2] )f»( 25, 18, 3.108) to (5, 8, 3.088) and another concentration around (19, 14, 3.074). 
In order to compute the volume density of the LAEs, we convolved the spatial distribution 
of the 56 LAEs with a Gaussian kernel with a — A Mpc, which is comparable to the redshift 



We note that there are differences between the redshifts of Lya emission- lines and those of other nebular 
lines which are expected to be better tracers of the systemic redshift of galaxies (Adelberger et al. 2003). 
Since neutral hydrogen and dust in galaxies tend to absorb the blue side of the Lya emission-lines, the peak 
of Lya emission- lines apparently shifts to higher redshifts. According to Adelberger et al. (2003), the excess 
and the rms scatter of redshifts for Lya emission-line for LAEs are 310 km s~^ and 250 km s~^, respectively. 
This scatter is smaller than the predicted peculiar velocity dispersion. 



- 5 - 



dispersion due to the predicted peculiar velocity dispersions, Av ~ 400 km s~^. By drawing 
the projected contour of volume density of 2.0 x 10~^ Mpc^'^, wc identified three filaments 
connecting with each other with the intersection at around (16, 11, 3.094). The length of 
each filament is about 30 Mpc and the width is about 10 Mpc in comoving scale. This is 
the largest coherent filamentary structure mapped in three dimensional space at 2; > 3. In 
the central 8.7' x 8.9' region, we also plot in Figure 1 the 21 LBGs in Steidel et al. (2003) 
whose Lya emission redshifts lie in the same redshift range of our narrow-band imaging 
observations, z — 3.054 — 3.120. These LBGs seem to be concentrated near the intersection 
of the filaments of LAEs. 

Although our spectroscopic observations are not complete yet, we tried to constrain the 
three-dimensional over-density of these filaments using the same Gaussian kernel above. The 
highest number density of the LAEs is ^ 6.0±2.4 x 10~^ Mpc~^. The average number density 
along the filaments is ^ 3.0 x 10~^ Mpc~^ while the average number density of 283 LAEs in 
the whole volume sampled by the narrow band filter, 1.4 x 10^ Mpc^, is 2.0 x 10~^ Mpc~^. 
Note, however, that the number density estimated from the spectroscopy is the lower limit 
because of the incompleteness in our redshift measurements. If we assume that the remaining 
66 LAEs, which are in the fields of 6 slit masks but not considered in the present analysis, 
have similar spatial distribution, the real number density of the LAEs in the filament would 
be higher by a factor of 2 and would be 3 times as large as the average value of the entire 
field. 

In Figure 2, wc show the composite spectrum of the 56 LAEs which arc shifted into the 
rest frame using their Lya emission redshifts. There is no evidence of [OIII] emission-lines 
for [Oil] emitters ai z ^ 0.33 in this spectrum. The rest frame EW of Lya emission-fine 
of the spectrum is about 60 A, which is roughly consistent with the EW estimated from 
our narrow-band imaging observations. The continuum fiux blueward of the Lya line is 
~20% dimmer than that redward of the line. However, this value is small compared with 
the continuum depression found in high S/N composite spectra of LBGs (e.g. Shapley et al. 
2003). This may be due to the poor quality of our spectra blueward of Lya emission-line 
because of the rapid decrease of the transmittance of our SY47 filter at 4600-4900A. The 
spectral resolution of 600 km s~^ is too large to resolve the line profile of Lya, since the typical 
FWHM of the Lya emission fine of LAEs at high redshifts is 300 km s~^ (e.g. Venemans et 
al. 2005). We could not find any significant evidence of other emission or absorption fines 
in the composite spectrum. The observed CIVA1549/Lya ratio is smaller than 0.1. This 
suggests that narrow-line AGNs do not dominate these LAEs; the CIVA1549/Lya ratio is 
typically 0.2 for narrow-line AGNs at 2; ~ 3 (Steidel et al. 2002). Thus the Lya emission of 
the LAEs is likely to originate from their star-formation activities. 



-6- 



A cosmological numerical simulation of galaxy formation including gas particles suggests 
that the average star-formation rate (SFR) of galaxies at z = 3 is nearly independent of 
their environment while a few galaxies with very large SFR may exist in the highest density 
regions (Keres et al. 2005). We estimated the SFR of the LAEs from their UV continuum 
luminosities. The left panel of Figure 3 shows the SFR of all the 283 LAEs as a function 
of projected surface number density, while the right panel shows the SFR as a function of 
volume density for the 56 LAEs with spectroscopic redshifts. We do not find any significant 
correlation. We checked that the dust reddening correction using their continuum colors does 
not change this trend. The environmental effect thus seems to be weak for the high-redshift 
LAEs, which is in fact consistent with the prediction of the simulation in Keres et al. (2005). 

There are lines of evidence that the LABI is a very massive galaxy with intensive star- 
formation activities. Bright submillimeter source was detected at the position of LABI 
and the SFR estimated from the submiUimeter fiux is extremely large (~ 1000 Mq yr~^) 
(Chapman et al. 2001). Not only the dust emission but the CO emission line was detected 
in LABI, which implies that the large amount of molecular gas also exists (Chapman et al. 
2004). The large mass of the host dark matter halo (~ 10^^ M©) was also suggested for LABI 
from the velocity dispersion and the physical extent of Lya emission-line, by assuming that 
the gas is bound within its gravitational potential (Bower et al. 2004). LAB2 looks similar 
to LABI on the Lya image, and it is fikely that the two giant LABs are very massive galaxies 
in their forming phase. Accordingly, it is interesting to see their location with respect to the 
filamentary structure. 

We measured the redshifts of LABI and LAB2 at their surface brightness peaks of 
the emission line. Their redshifts, z — 3.102 for LABI and z — 3.103 for LAB2, indicate 
that they are located near the intersection of the three filaments (Figure 1). Cosmological 
simulations predict that the intersections of large-scale filaments in the early Universe evolve 
into the present day massive clusters of galaxies. Thus, we can reasonably speculate that the 
two LABs may be progenitors of very massive galaxies near the center of a massive cluster. 
The smaller LABs of M04 are also concentrated near the position of the intersection in the 
projected distribution. It would be interesting to investigate by future observations whether 
or not the smaller LABs are preferentially located at the intersection of filaments in three 
dimensional space. 

Wc thank the anonymous referee for useful comments which have significantly improved 
the paper. We thank the staff of the Subaru Telescope for their assistance with our ob- 
servations. The research of T.Y. is partially supported by the grants-in-aid for scientific 
research of the Ministry of Education, Culture, Sports, Science, and Technology (14540234 
and 17540224). 
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Fig. 1. — {Bottom right panel) The sky map of the 283 candidate LAEs detected in Hayashino 
et al. (2004). The green hne shows the average local surface density of LAEs in this field 
(see the text). Cyan squares show two giant LABs. Blue circles show the field of view of 
6 masks. Blue, green, and red points show the LAEs at z = 3.05 — 3.08, z = 3.08 — 3.10, 
and ai z = 3.10 — 3.12, respectively. The triangles show the LBGs in the SSA22a field (the 
red box, 8.7' x 8.9', Steidel et al. 2003). {Top right and bottom left panels) The redshift 
space distribution of 56 LAEs with spectroscopic redshifts. The red line shows the projected 
contour of the local volume density of LAEs of 2 x 10~^ Mpc~^ (see the text). The predicted 
peculiar velocity dispersion of 400 km s~^ is shown by red arrows. The dotted lines show the 
redshift range sampled within > 50% of the peak transmittance of our narrow-band filter. 
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Fig. 2. — The median composite spectrum of the 56 LAEs shown in the rest frame. The 
vertical short dashed hnes show [OIII]AA4959,5007 for [OII]A3727 emitters at z = 0.33. The 
vertical dotted lines show wavelengths of strong interstellar absorption lines in Lyman break 
galaxies. The spectrum have been smoothed by a boxcar kernel of width 2.5A (the spectral 
resolution in the rest frame) . The long dashed hne shows the typical sky spectrum for LAEs 
at ^ = 3.09. 



- 11 - 



^ 100 
I 10 



> 



at 




0.1 



-1.5 -1 -0.5 

surface density logjoC-S / Mpc-2) 



100 



10 



> 

00 



0.1 



□ 

□ B- 



□ 
□ 

] 



(b) 



-3 -2.5 -2 

3d density logioCPsD / Mpc-^) 



Fig. 3. — The distribution of SFR estimated from UV continuum flux as a function of 
the projected local surface density S of LAEs {{a), left panel) and three-dimensional density 
PsD^ih) iTight panel). The squares show LAEs with spectroscopic redshifts and the points 
show LAEs without spectroscopic redshifts. The three dashed hues show the 75th percentile, 
the median, and the 25th percentile from the top to bottom. The three-dimensional density 
P3D is not corrected for the incompleteness in our redshift measurement. 



